Atmospheric 
Chemistry 
and Physics 


Atmos. Chem. Phys., 12, 5737-5753, 2012 
www. atmos-chem-phy s .net/ 1 2/5 7 37/20 1 2/ 
doi:10.5194/acp-12-5737-2012 
© Author(s) 2012. CC Attribution 3.0 License. 




Development of a climate record of tropospheric and stratospheric 
column ozone from satellite remote sensing: evidence of an early 
recovery of global stratospheric ozone 

J. R. Ziemke 1 ’ 2 and S. Chandra 3 

Goddard Earth and Sciences Technology and Research, Morgan State University, Baltimore, Maryland, USA 
2 NASA Goddard Space Flight Center, Greenbelt, Maryland, USA 

3 Goddard Earth Sciences and Technology, University of Maryland Baltimore County, Baltimore, Maryland, USA 
Correspondence to: J. R. Ziemke (jerald.r.ziemke@nasa.gov) 

Received: 17 October 201 1 - Published in Atmos. Chem. Phys. Discuss.: 30 January 2012 
Revised: 6 April 2012 - Accepted: 12 June 2012 - Published: 3 July 2012 


Abstract. Ozone data beginning October 2004 from the Aura 
Ozone Monitoring Instrument (OMI) and Aura Microwave 
Limb Sounder (MLS) are used to evaluate the accuracy of the 
Cloud Slicing technique in effort to develop long data records 
of tropospheric and stratospheric ozone and for studying their 
long-term changes. Using this technique, we have produced 
a 32-yr (1979-2010) long record of tropospheric and strato- 
spheric column ozone from the combined Total Ozone Map- 
ping Spectrometer (TOMS) and OMI. Analyses of these time 
series suggest that the quasi-biennial oscillation (QBO) is the 
dominant source of inter-annual variability of stratospheric 
ozone and is clearest in the Southern Hemisphere during the 
Aura time record with related inter-annual changes of 30- 
40 Dobson Units. Tropospheric ozone for the long record 
also indicates a QBO signal in the tropics with peak-to-peak 
changes varying from 2 to 7 DU. The most important result 
from our study is that global stratospheric ozone indicates 
signature of a recovery occurring with ozone abundance now 
approaching the levels of year 1980 and earlier. The nega- 
tive trends in stratospheric ozone in both hemispheres during 
the first 15 yr of the record are now positive over the last 
15 yr and with nearly equal magnitudes. This turnaround in 
stratospheric ozone loss is occurring about 20 yr earlier than 
predicted by many chemistry climate models. This suggests 
that the Montreal Protocol which was first signed in 1987 as 
an international agreement to reduce ozone destroying sub- 
stances is working well and perhaps better than anticipated. 


1 Introduction 

The measurement of tropospheric ozone from satellite re- 
mote sensing began with the studies by Fishman and 
Larsen (1987) and Fishman et al. (1990). In those studies 
tropospheric column ozone (TCO) was derived by subtract- 
ing co-located stratospheric column ozone (SCO) from to- 
tal column ozone. SCO was calculated from Stratospheric 
Aerosol and Gas Experiment (SAGE) ozone profile measure- 
ments while Total Ozone Mapping Spectrometer (TOMS) 
measurements were used for total column ozone. The method 
was named the tropospheric ozone residual (TOR) technique 
and was effective in producing global maps of tropospheric 
ozone albeit limited primarily to a seasonal climatology due 
to sparse data sampling for SAGE because of its occultation 
technique. Chandra et al. (2003) used the greater spatial and 
temporal coverage of stratospheric ozone from the Upper 
Atmosphere Research Satellite (UARS) Microwave Limb 
Sounder (MLS) in combination with TOMS total ozone to 
improve monthly maps of TOR. More recently Ziemke et 
al. (2006) and Schoeberl et al. (2007) used two separate 
residual algorithms combining Aura Ozone Monitoring In- 
strument (OMI) and Aura MLS measurements to produce 
global fields of TOR with improved temporal and spatial cov- 
erage beyond previous studies. 

An alternative residual technique is the convective-cloud 
differential (CCD) method (Ziemke et al., 1998) which uses 
total column ozone and SCO from the same satellite instru- 
ment. The CCD method is the most basic form of “cloud 
slicing” for measuring tropospheric and stratospheric ozone 
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(e.g., Ziemke et al., 2001, 2009, and references therein). 
Gridded measurements from the CCD method are restricted 
to the tropics; however mid-to-high latitude measurements 
may be obtained with this method and are given generally as 
zonal averages over the Pacific. An advantage with the CCD 
method is that by using a single instrument the potential er- 
rors arising from residual differencing are largely alleviated. 
This is important for deriving a long record of TCO and SCO 
from several combined instruments. 

The CCD technique takes advantage of the fact that UV- 
measuring instruments such as Total Ozone Mapping Spec- 
trometer (TOMS), Global Ozone Monitoring Experiment 
(GOME), and OMI do not measure ozone lying below op- 
tically thick clouds. The CCD method assumes that one can 
make an accurate estimate of SCO using highly reflecting 
deep convective clouds (reflectivity R > 0.8) in the tropical 
Pacific region. It was shown by Ziemke et al. (2009) from 
combined Aura OMI and MLS measurements that ozone 
concentrations inside thick clouds in the Pacific are small to 
within a few ppbv. The CCD method works well for deriving 
SCO in the Pacific because the above-cloud column ozone 
measured by OMI for thick clouds is essentially equivalent 
to SCO to within 1-2 DU. The CCD method of Ziemke et 
al. (1998) further assumes that SCO is invariant along longi- 
tude in the tropics. With these assumptions, TCO in tropical 
latitudes can be calculated at any given location by differenc- 
ing low reflectivity ( R <0.2) total column ozone and high 
reflectivity ( R > 0.8) SCO from the Pacific region within the 
same latitude range. Highly reflecting clouds generally do 
not have physical cloud heights reaching tropopause altitude 
and the column ozone above the cloud may vary consider- 
ably even when R > 0.8. As a practical solution, SCO in the 
Pacific is calculated using only minimum values of above- 
cloud column ozone in each gridded region. These minimum 
values are then averaged over the eastern and western tropi- 
cal Pacific to derive a single measurement of SCO. 

The validity of these assumptions has been established by 
comparing TCO and SCO derived from the CCD method 
with TCO and SCO derived from ozonesondes and satel- 
lite data including TOMS/MLS, TOMS/Halogen Occultation 
Experiment instrument (HALOE), and TOMS/SAGE com- 
binations (e.g., Ziemke et al., 1998, 2005; Chandra et al., 
2003). Application of the CCD method from TOMS data was 
extended over the Pacific region into the middle and high lat- 
itudes between 50° S and 60° N (Ziemke et al., 2005) for the 
time record 1979-2003. 

The objective of this study is to evaluate strengths and 
weaknesses of the CCD technique by incorporating Aura 
OMI and MLS ozone and to study inter-annual and decadal 
changes in tropospheric and stratospheric ozone. The Aura 
MLS measurements of SCO provide a stringent test of CCD 
SCO when compared to sparse satellite occultation measure- 
ments in previous studies. This is important for establishing 
a long record (1979-present) of CCD TCO and SCO from 
combined OMI and TOMS instrument measurements. Long 


record data products are important for addressing issues in- 
cluding trends in tropospheric ozone and pollution, strato- 
spheric ozone depletion, and decadal changes in the global 
Brewer Dobson Circulation. 

In the following, Sect. 2 discusses the Aura OMI and MLS 
measurements, Sect. 3 describes the CCD and OMI/MLS 
residual methodologies, Sect. 4 discusses inter-annual vari- 
ability of tropospheric and stratospheric ozone, Sect. 5 dis- 
cusses some results from the TOMS+OMI 1979-2010 ex- 
tended ozone dataset, and finally Sect. 6 provides a summary. 


2 Aura OMI and MLS ozone measurements 

The OMI and MLS are two out of a total of four instru- 
ments onboard the Aura spacecraft which is flown in a sun- 
synchronous polar orbit at 705 km altitude with a 98.2° in- 
clination. Aura was launched in July 2004 and has been pro- 
viding measurements since August 2004 to the present. The 
spacecraft has an equatorial crossing time of 01:45 p.m. (as- 
cending node) with around 98.8 min per orbit (14.6 orbits per 
day on average). Schoeberl et al. (2006) provide an overview 
of the EOS Aura mission and discuss the various measure- 
ments from the four Aura instruments. 

The OMI is a nadir-scanner which at visible (350-500 nm) 
and UV wavelength channels (UV-1: 270-3 14 nm; UV-2: 
306-380 nm) detects backscattered solar radiance to mea- 
sure column ozone with near global coverage (aside from 
polar night latitudes) over the Earth with a resolution 
of 13 x 24 km at absolute nadir. Aside from ozone, OMI 
also measures Optical Centroid Cloud Pressure (OCCP), 
aerosols, NO 2 , SO 2 , HCHO, and several other trace gases 
in the troposphere and stratosphere (Levelt et al., 2006). 
Measurements of ozone from OMI are determined using the 
OMT03 v8.5 algorithm which is an extension of the TOMS 
v8 algorithm. A description of the TOMS v8 algorithm may 
be obtained from the TOMS V8 CD DVD ROM, or from the 
OMI Algorithm Theoretical Basis Document (ATBD) from 
the TOMS web page (http://toms.gsfc.nasa.gov/version8/ 
v8toms_atbd.pdf). One difference between the TOMS v8 and 
the OMT03 v8.5 algorithms is the treatment of clouds. The 
TOMS v8 and earlier versions of OMT03 use a cloud pres- 
sure climatology based on thermal infrared cloud-top pres- 
sures, whereas OMT03 v8.5 uses in situ OCCP derived with 
OMI by the rotational Raman scattering method. 

SCO is calculated for the OMI/MLS residual method using 
MLS v2.2 ozone. The MLS instrument is a thermal-emission 
microwave limb sounder that measures vertical profiles of 
mesospheric, stratospheric, and upper tropospheric temper- 
ature, ozone, and several other constituents from limb scans 
taken in the direction ahead of the Aura satellite orbital track. 
The MLS profile measurements are made about 7 min before 
OMI views the same location during ascending (daytime) 
orbital tracks. These we refer to as “co-located” measure- 
ments between OMI and MLS. MLS also measures ozone 
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Tropospheric Ozone Measured 
by OMI/MLS and OMI CCD 



(CCD: "Convective Cloud Differential - i.e., most basic form of "Cloud Slicing") 

Fig. 1. Schematic illustration of the OMI/MLS tropospheric ozone 
residual method of Ziemke et al. (2006) (OMI total column ozone 
minus MLS stratospheric column ozone) and the convective-cloud 
differential (CCD) residual method of Ziemke et al. (1998) (OMI 
total column ozone minus OMI above-cloud column ozone). 


and other atmospheric constituents for descending nighttime 
orbits which on a given day can be up to ±12 h different in 
time from OMI daytime measurements. With combined as- 
cending and descending nodes MLS makes around 3500 ver- 
tical profile measurements over the Earth per day. This study 
includes only the ascending orbit co-located data from MLS 
for deriving SCO. Details regarding the instrument includ- 
ing spectrometers, spectral channels, calibration, and other 
topics are discussed by Waters et al. (2006) and in related 
papers in the same journal. Froidevaux et al. (2008, and per- 
sonal communication, 2011) provides validation results for 
MLS v2.2 measurements of ozone and other trace gases. At 
the present time an MLS v3.3 data product is also provided 
to the science community. While v2.2 retrieval has 37 pres- 
sure levels, v3.3 has 55 pressure levels and other improve- 
ments; however, v3.3 also has more outliers/missing data in 
the ozone measurements than with v2.2. Our analysis of SCO 
derived from MLS shows that there is little difference be- 
tween using v2.2 or v3.3 other than a small systematic off- 
set (v3.3 minus v2.2) of about 2.5 DU. Information regarding 
MLS v3.3 ozone measurements is available online from the 
NASA Data and Information Services Center (http://disc.sci. 
gsfc.nasa.gov/gesNews/mls_new_data_version_release). 


3 Retrieval methodologies and comparison of results 

Two residual methods are used for deriving TCO and 
SCO from Aura OMI and MLS measurements. One is the 
OMI/MLS approach of Ziemke et al. (2006) and the other 
is the CCD method of Ziemke et al. (1998). Figure 1 is a 
schematic diagram illustrating both of these residual tech- 


niques in the tropics where the tropopause is typically ~16- 
17 km altitude year-round. 

3.1 The OMI/MLS residual method 

For the OMI/MLS residual method in Fig. 1, SCO is de- 
rived from vertically integrated MLS ozone profiles which 
are subtracted from OMI near clear- sky (radiative cloud frac- 
tions < 30 %) total column ozone to derive TCO. Tropopause 
pressure, which separates tropospheric from stratospheric 
column ozone comes from National Centers for Environ- 
mental Prediction (NCEP) using the World Meteorological 
Organization (WMO) 2 Kkm -1 lapse rate tropopause defini- 
tion. SCO from MLS is determined by pressure integration 
of ozone volume mixing ratio profiles from 0.0215 hPa down 
to the NCEP tropopause. The MLS ozone profile measure- 
ments were linearly interpolated in log-pressure to the exist- 
ing NCEP tropopause pressure to derive SCO. MLS SCO (in 
Dobson Units, DU; 1 DU = 2.69 x 10 20 molecules m -2 ) was 
determined by standard log-pressure integration of ozone 
volume mixing ratio: SCO = 0.79 hPa^ ^P • d In P , 
where X is ozone volume mixing ratio in units ppbv and P is 
pressure in units hPa. The recommended range for scientific 
analysis of MLS v2.2 ozone profiles is 0.0215-215 hPa. As 
was done by Ziemke et al. (2006), nearly global SCO from 
MLS for each day was achieved by including ozone retrievals 
down to 316 hPa. 

MLS SCO data were binned to 1° latitude x 1.25° lon- 
gitude to be compatible with OMI level- 3 (L3) gridded to- 
tal column ozone. Tropopause pressures from NCEP anal- 
yses were re-binned to this same resolution from a coarser 
2.5 x 2.5° gridding. It is noted for MLS limb measurements 
that the horizontal optical path is about 300 km which is 
larger than the horizontal size of OMI L3 gridded data, but is 
comparable to the size of original NCEP gridded measure- 
ments. To derive a high density SCO field we have used 
the two-step spatial interpolation of Ziemke et al. (2006). 
The interpolation for SCO includes first a moving 2-D (lati- 
tude/longitude) Gaussian window along daytime orbit to fill 
in intermittent gaps along-track for MLS SCO, followed sec- 
ondly by a linear interpolation along longitude between ex- 
isting MLS SCO data. This interpolation approach preserves 
the along-track measurements of SCO from MLS at all lati- 
tudes. NCEP measurements of tropopause pressure were re- 
binned to the same 1° latitude x 1.25° longitude resolution. 
Following the determination of SCO and TCO at 1 x 1.25° 
resolution, the SCO data were averaged in 5 x 5° bins to be 
compatible with OMI SCO from the CCD method. 

3.2 The CCD residual method 

The CCD method illustrated in Fig. 1 uses a collective en- 
semble of high reflectivity (i.e., reflectivity R > 0.8) above- 
cloud column ozone measurements from OMI within speci- 
fied gridded regions (here 5° latitude x 5° longitude boxes) 
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to determine SCO. Cloud tops that lie lower in the tropo- 
sphere will have larger above-cloud column ozone amounts 
provided that spatial variability of ozone is small over the re- 
gion. High reflectivity scenes represent generally deep con- 
vective cloud systems, however their physical cloud tops may 
not generally reach tropopause altitude and column ozone 
above the cloud may vary considerably even with R > 0.8. 
As a practical approach, SCO in the Pacific with the original 
CCD algorithm of Ziemke et al. (1998) is estimated statisti- 
cally within each gridded region using only minimum values 
of above-cloud column ozone. This is done within each re- 
gion by subtracting 2 a from the mean value where a is calcu- 
lated RMS of the ensemble. All calculations are made each 
day with an absolute minimum number of R > 0.8 level-2 
footprint scenes in each 5x5° region chosen as 10. 

The minimum ensemble values are associated with deep 
convective clouds with small ozone concentration lying be- 
tween the OCCP and the tropopause. Evidence for the latter 
effect in the tropical Pacific was shown by Kley et al. (1996) 
during the Central Equatorial Pacific Experiment (CEPEX). 
Their study indicated from ozonesondes near-zero ozone 
concentrations in the upper troposphere with the passing of 
deep convective cloud systems. Possible mechanisms stated 
for the exceedingly low ozone measurements in the upper 
troposphere included chemical destruction of ozone by yet 
unidentified reactions and vertical transport via convective 
clouds of very low ozone concentrations from the low tro- 
posphere. Low ozonesonde concentration in the upper tro- 
posphere is also described in later studies by Folkins et 
al. (2002) and Solomon et al. (2005) as an indicator of deep 
convection and associated vertical injection of low ozone 
amounts upwards from the boundary layer/low troposphere. 
Vomel and Diaz (2010) in a more recent paper suggest that 
the near-zero ozone concentrations in the upper troposphere 
by Kley et al. (1996) were biased low because of uncertain- 
ties in ozonesonde cell currents. Following a reprocessing of 
the CEPEX ozonesondes, Vomel and Diaz (2010) indicated 
that lowest ozonesonde concentrations in the upper tropo- 
sphere in the Pacific may be more typically ~ lOppbv. As- 
suming a constant mixing ratio of lOppbv lying between 
the tropopause and OMI OCCP (~ 300-500 hPa for mini- 
mum above-cloud ozone scenes) this is equivalent to 1-3 DU 
in column amount. Even with such non-zero ozone concen- 
trations lying between the tropopause and OCCP for thick 
clouds, the estimate of SCO in the tropical Pacific from the 
CCD method should still be accurate to about 1-3 DU in ab- 
solute numbers. 

3.3 Comparison of results 

In the following, comparisons are made between tropo- 
spheric ozone derived from the two techniques followed by 
further discussion of the CCD method and ozone associated 
with deep convective clouds. 


3.3.1 Comparisons of CCD and OMI/MLS gridded 
tropical TCO measurements 

The primary measurement derived from the CCD method of 
Ziemke et al. (1998) is gridded TCO in low tropical latitudes. 
We have evaluated the CCD TCO data from OMI by compar- 
ing with coincident OMI/MLS residual measurements. Fig- 
ure 2 shows TCO from the CCD method (left panels) with 
TCO from OMI/MLS (right panels) for four consecutive Oc- 
tober months beginning in 2004. It is well known that Oc- 
tober (like September and November) coincide each year 
with large enhancements of tropospheric ozone in the tropi- 
cal south Atlantic. This zonal wave-one pattern in TCO in the 
tropics is caused largely by a combination of effects from the 
Walker circulation and photochemical sources which include 
biomass burning and lightning (e.g., Sauvage et al., 2007, and 
references therein). 

There is inter- annual variability present in Fig. 2 in the 
Atlantic region where year-to-year differences are ~ 10 DU. 
October 2004 and 2006 coincided with two El Nino events 
whereas October 2005 was non-ENSO (i.e., neither El Nino 
nor La Lina conditions) and October 2007 coincided with a 
strong La Nina event. Although ENSO events are predomi- 
nantly an Indian Ocean/Pacific Ocean phenomenon, it is pos- 
sible that ENSO may have contributed in some way to the 
inter-annual variability of ozone in the Atlantic. It is also pos- 
sible that the year-to-year changes in ozone in Fig. 2 are not 
related to ENSO. There is evidence in Fig. 2 of reduction 
of tropical tropospheric ozone in years 2004 and 2006 re- 
lated to the quasi-biennial oscillation (QBO) (discussed later 
in Sect. 5.1). The analysis of inter- annual variability of TCO 
in the Atlantic is beyond the scope and purpose of this study. 

Figure 3 shows a scatter plot comparing the two TCO 
products accumulated over the entire six years in the tropical 
domain 15° S-15° N (as shown in Fig. 2). The ensemble av- 
erages and RMS values listed in Fig. 3 are comparable for the 
two products with average values of about 27 DU and 6 DU, 
respectively. The calculated RMS of CCD minus OMI/MLS 
difference time series is much smaller at around 3 DU with 
an overall correlation between the two datasets of 0.88. The 
conclusion from Figs. 2 and 3 is that these two gridded prod- 
ucts are essentially similar with an average of about 3 DU 
statistical RMS difference and nearly zero offset. 

3.3.2 Tropospheric ozone associated with deep 
convective clouds 

A key uncertainty of the CCD method is the abundance of 
ozone lying inside deep convective clouds. Figure 4 shows 
seasonal contour plots of the difference of OMI CCD mini- 
mum above-cloud column ozone minus MLS SCO averaged 
over the six years. We refer to these residual differences in 
Fig. 4 as “cloud ozone” which represents tropospheric col- 
umn ozone lying between the tropopause and OCCP un- 
der conditions of deep convection (i.e., scene reflectivity 
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Fig. 2. Tropospheric column ozone (in Dobson Units) derived from the CCD residual method (left panels) and the OMI/MLS residual method 
(right panels) for four consecutive October months beginning in 2004. Blue to red colors represent smallest to largest values, respectively. 


exceeding 0.8). Some measurements in Fig. 4 are negative 
and exceed —5 DU in high latitudes, particularly during win- 
ter months. The negative column amounts in Fig. 4 are ob- 
viously not correct and are caused by yet unresolved offset 
differences between MLS SCO and OMI above-cloud col- 
umn ozone in wintertime high latitudes where solar zenith 
angles are highest, often exceeding 75°. 

It is possible that column differences in Fig. 4 could in part 
involve ozone measurement errors over thick clouds from 
OMI. Liu et al. (2004) discusses ozone retrieval errors (posi- 
tive offset errors) in column ozone measurements over thick 
clouds caused by the assumption of opaque Lambertian cloud 
surfaces in the TOMS algorithm. Liu et al. (2004) invoked 
a radiative transfer model to investigate an older version 7 
TOMS ozone algorithm which used infrared cloud pressure 
climatology. In comparison, OMI ozone measurements in 
Fig. 4 use in situ OCCPs which vary with each footprint mea- 
surement and generally lie several hundred hPa lower in alti- 
tude than IR cloud pressures. It is difficult to assess potential 
OMI cloud ozone errors with the Liu et al. (2004) TOMS 
version 7 results. However, Ziemke et al. (2009) did show 
for tropical latitudes that the OMI ozone retrievals over thick 


clouds yielded consistent and well-calibrated ozone concen- 
trations and column amounts. 

Ideally, wherever the CCD method correctly derives local 
measurements of SCO the differences between the two ozone 
datasets in Fig. 4 should be zero. The differences in Fig. 4 
are generally small in the Pacific extending from the trop- 
ics to high latitudes, but there are offset differences as large 
as 10 DU or even 20 DU in the extra-tropics of both hemi- 
spheres in regions away from the Pacific. It will be shown 
that these offset differences are mostly recurring annual cycle 
features and that tropospheric and stratospheric ozone from 
the CCD method are more accurate after these measurements 
are deseasonalized. The small differences in the Pacific sec- 
tor motivate the use of the CCD technique in analyses to 
follow. 

The column ozone differences in Fig. 4 in the tropics 
are consistent with the tropical cloud-ozone concentrations 
measured by Ziemke et al. (2009); that is, low ozone in 
the western-eastern Pacific year-round with high ozone over 
northern and southern Africa for July and October months 
(e.g., their Figure 9). The method used to derive cloud ozone 
by Ziemke et al. (2009) was different than the method used 
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Fig. 3. Scatter plot of CCD versus OMI/MLS monthly mean gridded 
tropospheric column ozone (in Dobson Units) accumulated over the 
six-year record in the tropics. The measurements for both products 
are gridded at 5 x 5° latitude-longitude resolution and extend from 
central latitudes 12.5° S to 12.5° N along all longitudes. 


in Fig. 4. Ziemke et al. (2009) determined OMI above-cloud 
column ozone by using OMI in situ OCCP measurements. 
Above-cloud column ozone in Fig. 4 is instead determined 
from the CCD method which does not require cloud pres- 
sures. 

Ziemke et al. (2009) surmised that deep convection in 
the tropics could explain much of the measured low-to-high 
ozone concentrations associated with thick clouds. That is, 
low tropospheric ozone values in the Pacific and high tro- 
pospheric ozone values over Africa and South America for 
thick cloud scenes were hypothesized as being largely a man- 
ifestation of boundary layer ozone injected into the clouds 
from below. Ziemke et al. (2009) however could not dis- 
tinguish between boundary layer ozone injected into the 
thick clouds and background ambient ozone lying above 
OCCP as a mixture from other sources such as stratosphere- 
troposphere exchange (STE), lightning, and biomass burn- 
ing. Despite smaller influence from STE in the tropics it was 
not possible to distinguish cloud-injection from other sources 
for the measured ozone concentrations over deep convective 
clouds. 

In Fig. 4 it is not possible to assess the sources for the 
large positive cloud ozone columns for any region, espe- 
cially in the extra- tropics in and near tropospheric wind jets 
where there is large day-to-day spatial variability and uncer- 
tainties in column ozone amounts. Furthermore, the largest 
CCD minus MLS column ozone differences in Fig. 4 are 
not necessarily a direct result of local deep convection and 
large boundary layer ozone concentration, but instead could 
be caused by substantial amount of ozone lying above OCCP 


from other extra-tropical sources (such as STE) under condi- 
tions of low cloud height and high reflectivity. 

The OMI measurements in Fig. 4 include optically thick 
clouds with R > 0.8, but even optically thick clouds in many 
regions are not deeply convective and have geometrical cloud 
tops often lying well below tropopause. This condition ap- 
plies to most of the oceanic regions aside from the Inter- 
Tropical Convergence Zone (ITCZ). The large positive col- 
umn ozone amounts in the subtropics of both hemispheres in 
Fig. 4 coincide most generally with these regions were OCCP 
lies well below tropopause pressure for thick cloud scenes. 

It is possible that the large positive ozone differences in 
Fig. 4 can be explained as substantial ozone amounts ly- 
ing above thick cloud coming from various sources (such 
as STE) simply under the conditions of low cloud height 
and high reflectivity. One might predict that largest positive 
column ozone amounts in Fig. 4 will coincide with largest 
differences between tropopause and OCCP for thick cloud 
scenes, but it turns out that this is not the case. We have cal- 
culated OCCP minus tropopause pressure for R > 0.8 scenes 
and find that there is no significant positive correlation be- 
tween the pressure differences and the column ozone differ- 
ences in Figure 4. The pressure thickness comparisons do not 
contribute to understanding the ozone differences in Fig. 4 
for deep convective clouds and are not shown here. One 
would need to know the amount of surface/boundary layer 
ozone injected into thick clouds from below relative to other 
sources of ozone lying within and above the clouds. Under- 
standing the sources and magnitude of cloud ozone in Fig. 4 
would require extensive cloud/convection modeling which is 
beyond the scope of our study. 

Comparisons between ambient ozone and cloud ozone 

Line plots of 12-month annual cycles of CCD cloud ozone 
in the Pacific in 20° latitude bands from 60° S to 60° N are 
shown as solid curves in Fig. 5. Also plotted as dotted curves 
in Fig. 5 are corresponding TCO time series from OMI/MLS. 
We refer to these measurements as background “ambient” 
tropospheric ozone. Included in Fig. 5 for all time series are 
calculated standard error of mean numbers which provide a 
measurement proportional to inter-annual variability; as ex- 
ample, inter-annual variability for cloud ozone in high lati- 
tudes in Fig. 5 is around 3-4 time larger than in the tropics. 

The annual cycles and annual mean values for cloud ozone 
measurements in the tropics in Fig. 5 are small at only about 
1-3 DU, however annual cycles and annual means can ex- 
ceed 5 DU in the extra- tropics. Ambient ozone in Fig. 5 
maximizes in spring-summer months (March-July) in the 
Northern Hemisphere and the spring months (September- 
November) in the Southern Hemisphere. The cloud ozone in 
Fig. 5 is always substantially smaller than ambient ozone at 
all latitude ranges, yet within the RMS uncertainties the an- 
nual cycles are generally similar for both column amounts. 
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Fig. 4. Three-month seasonal averages (December-February, March-May, June-August, and September-November) of CCD above-cloud 
column ozone minus MLS stratospheric column ozone over the time record 2004-2010. Contour numbers represent Dobson Units. The 
colors violet/blue to orange represent negative to positive values, respectively. We denote these measurements in this figure as “cloud ozone”, 
the amount of column ozone lying between the tropopause and cloud OCCP under conditions of deep convection. White regions denote areas 
where there is insufficient number of deep convective clouds for the CCD method. 
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Fig. 5. Solid curves: line plots of CCD 12-month annual cycles of cloud ozone averaged within six non-overlapping 20° latitude bands in the 
Pacific. The latitude bands are (from upper left to lower right): 40-60° N, 20-40° N, 0-20° N, 0-20° S, 20-40° S, and 40-60° S. The Pacific 
averaging is for the longitude domain 120° W-120° E about the dateline. Dotted curves: same as solid curves but for background ambient 
tropospheric column ozone from OMI/MLS residual. These annual cycle time series are all derived by averaging together data for similar 
months over the six years. The vertical bars represent calculated ±1 o RMS standard error of mean. 
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Fig. 6. Similar to Fig. 5 but for zonal mean measurements. 


West Coast USA 25N-45N,120W-150W 


Mediterranean 20N-40N,0-30E 
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Fig. 7. Similar to Figs. 5 and 6 but instead for six broad regions of the globe where measured CCD cloud ozone is large with peak abun- 
dances of ~ 15-20 DU. These six regions are from upper left to lower right (compare these regions with Fig. 4): West Coast of the USA, 
Mediterranean, Southeast Asia, West Coast of South America, Southern Africa, and West Coast of Australia. 


Figure 6 is the same as Fig. 5 except that the time series 
were derived for zonal means rather than Pacific means. The 
conclusions for zonal means are similar to those discussed 
for Pacific means except that cloud ozone annual means are 
larger at most latitudes for zonal mean measurements. Fig- 


ure 7 shows more annual cycle line plot comparisons, but in- 
stead for six extra-tropical regions where cloud ozone often 
exceeds 5-20 DU (i.e., exceeding 50% of ambient ozone in 
some months). It is inferred that these extra- tropical regions 
may be more polluted with boundary layer ozone injected 
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OMI CCD Pac SCO (DU) Monthly Means 



MLS Pac SCO (DU) Monthly Means 



Fig. 8. Latitude versus month contour diagrams of monthly mean 
stratospheric column ozone (in Dobson Units) averaged over the 
Pacific (120° W-120° E) from OMI CCD (top) and MLS (bottom). 
Dark to light shading designates smaller to higher column amounts, 
respectively. The contour values for both diagrams begin at 200 
Dobson units and increment by 20 Dobson Units. 

into the thick clouds from below and are more likely to have 
higher ozone concentrations in the upper troposphere in and 
above clouds coming from sources including STE and pollu- 
tion. The annual cycles for cloud ozone and ambient ozone 
are not correlated for Southeast Asia and southern Africa 
in Fig. 7. For southern Africa, ambient ozone maximizes in 
September-October (same as the other two Southern Hemi- 
sphere regions in the figure), whereas cloud ozone maxi- 
mizes earlier around August; August is a peak month for 
biomass burning in the southern Africa region. 

The conclusions from Figs. 4-7 are that Pacific means and 
zonal means have nearly the same geophysical signatures in 
annual cycles for both the cloud ozone and ambient ozone, 
and also that most of the offset differences in Fig. 4 are re- 
curring annual cycle features. Later in Sect. 4 we will show 
that after removing annual cycles from the data that inter- 
annual variability for the Pacific mean and zonal mean data 
products of both TCO and SCO from the CCD method are 
within a few DU from OMI/MLS residual ozone, not just in 
the tropics but extending to high latitudes. 

4 Inter-annual variability of tropospheric and 
stratospheric ozone 

In the study by Ziemke et al. (2005) the CCD measurements 
of TCO and SCO from combined Nimbus 7 and Earth Probe 
TOMS were extended to middle and high latitudes in the Pa- 


MLS Pac SCO (DU) Deseas Monthly Means 



Month 


OMI CCD Pac SCO (DU) Deseas Monthly Means 



Month 

Fig. 9. Contour diagrams of the same stratospheric column ozone 
(in Dobson Units) plotted in Fig. 8, but instead with all of the data 
deseasonalized. Solid contours represent positive ozone anoma- 
lies and begin with +2 Dobson Units increasing in increments of 
+2 Dobson Units. Dotted contours represent negative ozone anoma- 
lies and begin with —2 Dobson Units decreasing in increments of 
—2 Dobson Units. 


cific under the assumption that deep convective clouds with 
low boundary layer ozone conditions persist outside the trop- 
ics just as they do in the tropics. The CCD measurements 
of SCO were tested against SAGE II SCO for the time pe- 
riod 1984-2003. Despite sparse SAGE measurements those 
comparisons showed that SAGE and CCD SCO in the extra- 
tropics compared well in annual means and even better (to 
2-4 DU differences) from the tropics to high latitudes when 
comparisons were made for summer months only. 

We use the Aura MLS measurements of SCO which have 
greater spatial and temporal coverage than SAGE to evalu- 
ate how well the CCD method works outside tropical lati- 
tudes. Figure 8 compares latitude versus month contour di- 
agrams of SCO averaged over the Pacific (120° W-120° E) 
from OMI CCD (top) and MLS (bottom). SCO in middle 
and high latitudes for either OMI or MLS in Fig. 8 is largest 
in both hemispheres from winter into spring with large 80- 
100 DU peak-to-peak annual-cycle changes. There is indica- 
tion of large inter-annual variability in both data sets from 
the tropics to high latitudes which includes an apparent QBO 
signal in the tropics with a period of about two years. 

The SCO data in Fig. 8 were further deseasonalized to 
evaluate inter- annual changes (Fig. 9). Deseasonalization 
was accomplished by subtracting for each month a global 
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monthly mean climatology value (determined by averaging 
similar months over the six-year record). The dominant inter- 
annual signature in Fig. 9 is the QBO which during the Aura 
record has about a 24-month period and is characteristically 
centered about the equator extending to high latitudes of both 
hemispheres. 

The tropically driven QBO induces a global secondary cir- 
culation with opposite vertical wind fields between the trop- 
ics and extra-tropics (e.g., Andrews et al., 1987). Subsidence 
(ascent) of stratospheric air mass in equatorial latitudes as- 
sociated with the QBO-induced circulation coincides with 
ascent (subsidence) of stratospheric air mass in the extra- 
tropics. In the tropics the subsidence of air mass driven by the 
QBO during the descending westerly phase (i.e., descend- 
ing eastward zonal winds in the low-mid stratosphere) is 
seen in Fig. 9 as anomalous increases in tropical SCO that 
maximize around the months May-October for years 2006 
and 2008. The QBO-induced down- welling circulation in the 
tropics coincides with opposite upwelling in the extra-tropics 
which is seen in Fig. 9 as anomalous reductions in SCO of 
10-20 DU in high latitudes of both hemispheres in winter- 
spring months. The extra-tropical QBO variability appears 
clearer in the Southern Hemisphere because of a synchronous 
phase coupling between the Brewer Dobson Circulation an- 
nual component and the QBO which during the Aura record 
had about a 24-month cycle. The synchronous phase cou- 
pling between annual cycle and QBO in ozone and other 
stratospheric trace gases was studied by Dunkerton (2001). 
Dunkerton (2001) employed HALOE satellite measurements 
for an earlier time record of similar length (i.e., 1992-1999) 
and cites several studies relating to synchronization between 
QBO and annual cycle variability. 

Toward the end of the record in Fig. 9 there is anomalously 
large SCO in the northern higher latitudes beginning around 
January 2010 for both OMI and MLS. These increases co- 
incide with ozone decreases in the tropics associated with 
the descending easterly phase of the QBO. A recent study 
by Steinbrecht et al. (2011) combined ozonesondes from 
Hohenpeissenberg (48° N, 11°E) and SCIAMACHY total 
ozone to conclude that these high values of ozone in 2010 
were among the largest on record in northern latitudes dur- 
ing the last 20-25 yr. Their study attributes these large ozone 
enhancements to a coupling between the QBO and the Arctic 
Oscillation and North Atlantic Oscillation with the latter two 
oscillations being in an unusually persistent negative phase. 

Time series of SCO corresponding to the data in Figs. 8 
and 9 are plotted in Fig. 10 in 10° latitude bands (indi- 
cated) for monthly means (left panels) and deseasonalized 
monthly means (right panels). Shown in all panels in Fig. 10 
are Pacific averages for CCD (solid curves) and MLS (dotted 
curves). Despite seasonally varying offsets up to 5-10 DU 
between CCD and MLS in the left panels in Fig. 10, when 
the data are deseasonalized the month-to-month differences 
between the two ozone datasets is reduced for all latitude 
ranges. The deseasonalized time series in Fig. 10 for both 


residual methods show that the QBO-related signal in the 
Southern Hemisphere high latitudes is associated with inter- 
annual changes of ~ 30 DU or greater. 

Pacific mean tropospheric ozone derived from the two 
residual methods is shown in Fig. 11. Figure 11 is the same 
as Fig. 10 but with TCO plotted rather than SCO. Shown in 
Fig. 1 1 are the original time series (left panels) and desea- 
sonalized time series (right panels) for the same five latitude 
bands. There are obvious offset differences between the two 
methods in the left panels in Fig. 1 1 which are up to ~ 10 DU 
in some months. The offset differences for TCO in Fig. 11 
are equal in magnitude (but of opposite sign) to the SCO off- 
set differences in Fig. 10 since the same OMI total column 
ozone is used for both TCO time series. The deseasonalized 
TCO time series in the right panels do not have these offsets 
and track each other reasonably well. 

Figures 9-1 1 suggest that inter-annual variability of Pa- 
cific mean SCO and TCO from the CCD method from the 
tropics to high latitudes compares closely to within a few DU 
with corresponding measurements from the OMI/MLS resid- 
ual method. We will show that zonally averaged CCD above- 
cloud column ozone also compares close to zonal mean SCO 
from MLS. 

Figure 12 shows contour diagrams of deseasonalized zonal 
mean SCO from MLS (top panel) and OMI CCD above- 
cloud column ozone (bottom panel). Comparison with Fig. 9 
suggests that SCO from the two methods are closer for zonal 
means than for Pacific means. 

Figure 13 plots temporal RMS values of the difference be- 
tween CCD and MLS deseasonalized SCO time series for 
Pacific means (asterisks) and zonal means (triangles) as a 
function of latitude. (The RMS amplitudes for Pacific means 
and zonal means in Fig. 13 were calculated from OMI mi- 
nus MLS differences of the data plotted in Figs. 9 and 12, 
respectively.) The RMS values in Fig. 13 for most latitudes 
are generally about 1-2 DU and up to ~ 2-3 DU at higher 
latitudes. For zonal means RMS differences are about 0.5- 
1 DU in the tropics to ~ 1.5 DU at mid-to-high latitudes. The 
RMS values in Fig. 13 for SCO are equivalent to the RMS 
amplitudes calculated for TCO. Zonal mean measurements 
of CCD deseasonalized SCO in Fig. 13 have slightly smaller 
RMS amplitudes than the Pacific mean measurements. The 
smaller RMS amplitudes for zonal mean measurements re- 
flects reduced noise because of more gridded measurements 
and also that the offsets present between CCD and MLS SCO 
in all regions including the Atlantic are largely recurring sea- 
sonal cycle variations which are largely eliminated following 
deseasonalization of the data. 

These comparisons between Aura OMI and MLS column 
ozone have given us greater confidence in the CCD measure- 
ments of tropospheric and stratospheric column ozone de- 
rived from previous TOMS retrievals. In the next section we 
discuss a long record of stratospheric and tropospheric ozone 
determined from combined TOMS and OMI records begin- 
ning 1979. 
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Fig. 10. (Left panels) Monthly averaged Pacific mean measurements of CCD SCO from OMI (solid curve) and MLS SCO (dotted curve) 
averaged over five indicated 10° latitude bands. The Pacific mean represents data averaged over the combined eastern and western Pacific 
(i.e., longitude range 120° W to 120° E about the dateline). All measurements are in Dobson Units. (Right panels) Same as left panels but 
with the data deseasonalized. 

5 The TOMS+OMI ozone dataset spheric ozone. In the near future we plan to combine the 

TOMS and OMI stratospheric and tropospheric ozone data 
We have developed an extended record of tropospheric and using only the single version 9 processed retrievals. We pro- 

stratospheric column ozone spanning 1979-2010 by combin- vide only a brief discussion of ozone trends and inter-annual 

ing TOMS and OMI CCD measurements. The OMI data are (QBO) variability in tropospheric ozone for the 1979-2010 

appended with the earlier TOMS data used by Ziemke et data record, 

al. (2005). It was noted in the data description section that 

the TOMS ozone uses the version 8 algorithm while OMI 5.1 Quasi-Biennial Oscillation signals in tropical ozone 
ozone uses the version 8.5 algorithm. Although there may be 

retrieval offsets existing between TOMS and OMI measure- Although it is well known that there exists a large QBO 

ments because of different algorithms for the separate instru- variability in stratospheric ozone, there is also evidence of 

ments, our analyses suggest that these offsets are not large a QBO in tropospheric ozone. Ziemke and Chandra (1999) 

and cannot be more than a few DU at most at any latitude. first detected a possible QBO in tropospheric ozone using 

Offsets of only a few DU will adversely affect evaluation of Nimbus-7 TOMS measurements for the 1979-1993 record, 

inter-annual variability and calculated trends in tropospheric A later study by Chandra et al. (2002) combined Nimbus- 

ozone, but such offsets have less relative impact for strato- 7 TOMS with Earth Probe TOMS and found a consistent 
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Fig. 11. (Left panels) Monthly averaged Pacific mean measurements of CCD TCO from OMI (solid curves) and OMI/MLS TCO (dotted 
curves) averaged over five indicated 10° latitude bands. The Pacific mean represents data averaged over the combined eastern and western 
Pacific (i.e., longitude range 120° W to 120° E about the dateline). All measurements are in Dobson Units. (Right panels) Same as left panels 
but with the data deaseasonalized. 


QBO in tropospheric ozone over a longer record of 1979- 
2000. The QBO in tropospheric ozone for these studies was 
clearer to detect in the Atlantic region as opposed to the 
Pacific as there is additional inter-annual variability in the 
Pacific related to ENSO events. The QBO in tropospheric 
ozone was shown to be negatively correlated with the QBO 
in SCO. It was postulated by Ziemke and Chandra (1999) 
that a possible source for the QBO signal could be an upper 
tropospheric ozone photochemistry response from the QBO 
in stratospheric ozone. However, the measured tropospheric 
signal was larger than predicted by photochemical model- 
ing and it was concluded by Ziemke and Chandra (1999) 
that the source for the QBO signal in tropospheric ozone 
was most likely of dynamical origin. Lee et al. (2010) have 
reached these same conclusions based upon balloon sonde 
data from the Southern Hemisphere Additional Ozonesondes 
(SHADOZ) network. Their analyses of the sonde tempera- 


ture and ozone measurements suggest that the QBO signal in 
tropospheric ozone extends down to about mid-troposphere 
and is of dynamical origin. 

Figure 14 plots 50hPa monthly zonal winds from Singa- 
pore (1° N, 104° E) (dotted curve) and deseasonalized tropo- 
spheric ozone (solid curve) averaged in the equatorial At- 
lantic (5° S-5° N, 60° W-60° E). There are two data gaps 
present in the time series in Fig. 14. The first gap is several 
years of non-existing measurements between Nimbus 7 and 
Earth Probe TOMS periods. The second gap in Fig. 14 repre- 
sents Earth Probe CCD measurements which have been con- 
servatively flagged as missing after year 2000 for question- 
able data quality. Tropospheric ozone in Fig. 14 was addi- 
tionally smoothed using a low-pass digital filter as described 
in the figure caption. 

Comparison of the 50hPa zonal winds and tropospheric 
ozone in Fig. 14 indicates a persistent negative correlative 


Atmos. Chem. Phys., 12, 5737-5753, 2012 


www.atmos-chem-phys.net/12/5737/2012/ 










J. R. Ziemke and S. Chandra: Development of a climate record 


5749 


MLS ZM SCO (DU) Deseas Monthly Means 



Month 


OMI CCD ZM SCO (DU) Deseas Monthly Means 



Month 

Fig. 12. Same as Fig. 9, but instead for zonal means rather than 
Pacific means. 


Temporal RMS Calculations 



Fig. 13. Calculated temporal RMS values (in Dobson Units) of the 
difference between CCD and MLS deseasonalized SCO time series 
for Pacific means (asterisks) and zonal means (triangles) as a func- 
tion of latitude. These RMS amplitudes were calculated from the 
difference between the data plotted in Figs. 9 and 12, respectively. 


relationship over the long 32-yr multi-instrument record. 
The QBO signal in tropospheric ozone for the 1979-1993 
Nimbus-7 TOMS period in Fig. 14 is larger than during the 
latter EP TOMS and OMI periods. Peak-to-peak differences 
in ozone for the Nimbus-7 record are ~ 4-7 DU compared 
to ~ 2-4 DU for the latter years. There are some years in 
Fig. 14 where the wind/ozone negative correlation relation is 
not evident. One case occurs around 1990-1991 (coinciding 



Fig. 14. Tropospheric column ozone from the CCD method in Dob- 
son Units (solid curve) plotted versus 50 hPa zonal winds from Sin- 
gapore (1° N, 104° E) in units ms -1 but divided by 10 for scaling 
(dotted curve). The ozone time series was averaged over the equa- 
torial Altantic region (5° S-5° N, 60° W-60° E). The ozone time 
series was deseasonalized and detrended, and then smoothed us- 
ing a recursive low-pass digital filter with one-half filter response at 
12-month period (filter response of about 0. 8-0.9 for 24-36 month 
QBO time periods). 


with the July 1991 Mt. Pinatubo eruption), and another case 
is in 1997-1998 during an intense tropical El Nino event. 

We conclude from the TOMS+OMI 32-yr record that there 
appears to be a persistent QBO signal in tropospheric ozone 
with peak-to-peak amplitudes varying from about 2 DU up 
to 7 DU. This evidence further strengthens claims from pre- 
vious studies of a QBO in tropospheric ozone. We note that 
Lee et al. (2010) from ozonesonde analyses found maximum 
QBO signal in tropospheric ozone profile measurements of 
about 8 ppbv (equivalent to ~2-3 DU in upper troposphere) 
which is not inconsistent with the signal amplitudes that we 
find from the TOMS+OMI combined record. 

5.2 Ozone trends 

Trends in tropospheric and stratospheric ozone were calcu- 
lated by Ziemke et al. (2005) for a 25-yr record (1979- 
2003) of TOMS CCD measurements in the Pacific. Ziemke 
et al. (2005) found that trends in tropospheric ozone were 
near zero most everywhere but marginally positive (~ 
+2 DU to +3 DU decade -1 ) in mid-latitudes of both hemi- 
spheres. Trends in stratospheric column ozone were also 
near zero in the tropics but large and negative (~ — 10 to 
— 14 DU decade -1 ) in the mid-high latitudes of both hemi- 
spheres. We have made similar calculations of trends us- 
ing the extended TOMS+OMI 32-yr record. The results pre- 
sented are limited to stratospheric ozone trends; tropospheric 
ozone trends for the 32-yr record are not substantially differ- 
ent from Ziemke et al. (2005) based on a shorter data record. 

Figure 15 shows calculated trends (both DU-decade -1 and 
percent-decade -1 units) in Pacific mean stratospheric col- 
umn ozone for the Nimbus-7 15-yr (1979-1993) beginning 
record (bottom dotted curve), and 15-yr (1996-2010) end- 
ing record (top dashed curve). For comparison, the aver- 
age of these two trend curves is also shown (middle solid 
curve). The trends in Fig. 15 for the two 15-yr datasets were 
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Fig. 15. Trends of CCD stratospheric column ozone for the 1979- 
1993 beginning 15-yr record (bottom dotted curve) and 1996-2010 
ending 15-yr record (top dashed curve). The trend units for the left 
vertical axis are DU decade - 1 with the right vertical axis units given 
as percent decade -1 . As comparison, trends for these two 15-yr 
records are averaged together (middle solid curve). Trends were de- 
rived for the two 15-yr records using a linear multivariae regression 
model (see text). All data were averaged over the Pacific (120° W 
to 120° E about dateline) within 5° latitude bands. Vertical bars in- 
dicate ±2(7 annual mean trend uncertainties calculated from the re- 
gression model. The trend uncertainty bars include an additional 
2 DU (2cr ) long-term inter-instrument calibration/drift uncertainty. 
Values of 2cr for the averaged trends were derived by taking the 
square root of the sum of the variances for the two 15-yr record 
variances x2. 


Column Ozone 60 S- 60 N Column Ozone 60 S- 60 N 




Fig. 16. (Left) Annual mean time series of MOD total column ozone 
(thick solid curve along top), TOMS+OMI CCD measurements of 
stratospheric column ozone (thick dotted curve along top), and so- 
lar F10.7 cm radio flux (thin solid curve along bottom) which has 
been rescaled for plotting. Stratospheric column ozone has had 31 
Dobson Units added for visualization in plotting with total column 
ozone. Both ozone time series represent area- weighted measure- 
ments lying between latitudes 60° S and 60° N. The MOD total 
column ozone time series was derived from zonal means whereas 
the CCD stratospheric column ozone time series was derived from 
Pacific means. (Right) Similar to ozone plotted in the left panel 
except that time averages have been removed and the 11 -yr solar- 
cycle signals in the two ozone time series have been extracted by 
linear regression. The regression model used to remove 11 -yr so- 
lar cycle variability in the two ozone time series is ( t is year in- 
dex): Ozone(t) = C - Solar (t) + s(t) where C and is a constant, So- 
lar is solar FI 0.7 cm annual mean time series (with time average 
removed) and e = model residual error. (The two derived residual 
series s(t) are plotted in the right panel.) A 3-yr running average 
was applied to both time series in the right panel for smoothing. 


calculated for monthly means using the multivariate regres- 
sion model of Ziemke et al. (2005). This model included re- 
gression fits for combined linear trend, seasonal cycle, QBO, 
solar cycle, and ENSO. The latitude range in Fig. 15 coin- 
cides with the trend figures shown by Ziemke et al. (2005) 
in which there are too few Nimbus 7 CCD monthly measure- 
ments pole- ward of 50° S for seasonal trend analysis. The 
two 15-yr beginning and ending periods were chosen based 
on plotting of time series such as in Fig. 16 (discussed later) 
which has a characteristic “V”-shaped time dependence after 
removing solar signal. 

It is evident in Fig. 15 that the strong negative trends 
in stratospheric ozone outside the tropics for the 15-yr be- 
ginning period (bottom dotted curve) are nearly negated by 
trends for the 15-yr ending record (top dashed curve). Most 
of the change in trends and global increase in ozone over the 
last 15yr lies in the Northern Hemisphere extending from 
the subtropics to high latitudes. Area weighting of the extra- 
tropical trend differences (i.e., ending minus beginning 15- 
yr record trends) in Fig. 15 for latitudes 50° S to 50° N in- 
dicates 2/3 of the change in stratospheric ozone lies in the 
Northern Hemisphere. Trends of 10 to 20 DU decade -1 in the 
mid-high latitudes in Fig. 15 correspond to trends of approxi- 


mately 3 to 6 % decade -1 , respectively. Inter-instrument cal- 
ibration offsets and/or drifts cannot explain the large changes 
in trends in the extra- tropics in Fig. 15 between the two time 
records. That is, calibration errors are not generally depen- 
dent of latitude, and although large trend changes are mea- 
sured in the extra-tropics in Fig. 15, trends in the tropics re- 
main near zero for both time records. 

A complicating factor in assessing global stratospheric 
ozone recovery and turnaround period of stratospheric ozone 
loss is solar cycle variability. Figure 16 (left panel) shows 
annual mean time series of CCD stratospheric ozone (stars) 
plotted with merged total column ozone from the NASA 
merged ozone datasets (MOD) webpage http://acdb-ext.gsfc. 
nasa.gov/Data_services/merged/. Both time series represent 
column ozone area averaged between 60° S and 60° N. The 
MOD total column ozone time series in Fig. 16 was de- 
rived from zonal means whereas the CCD stratospheric col- 
umn ozone time series was derived from Pacific means. It 
was shown Ziemke et al. (2005) that trends in total col- 
umn ozone and stratospheric column ozone are nearly iden- 
tical for zonal means and Pacific means either globally or as 
function of latitude. For CCD stratospheric ozone a constant 
3 1 DU was added to the time series for plotting with total 
column ozone. This 3 1 DU accounts for the abundance of 
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tropospheric ozone. Along the bottom in Fig. 16 as a proxy 
of solar UV variability is solar F 10.7 cm radio flux time se- 
ries (e.g., Jackman et al., 1996, and references therein). 

One cannot readily identify a turnaround period in the 
ozone time records in the left panel of Fig. 16 because of a 
dominant solar cycle present. In the right panel of Fig. 16 we 
have removed most solar cycle variability in both ozone time 
series using regression (discussed in figure caption). With so- 
lar variability removed we identify a turnaround period oc- 
curring in the mid- 1990’s for stratospheric and total column 
ozone. We can also conclude from the right panel in Fig. 16 
that tropospheric ozone does not indicate substantial decadal 
changes. 

Several previous studies have done similar analyses as in 
Figs. 15-16 but were based upon shorter ozone records. We 
mention some of these studies which have relevant impact 
for our results. Weatherhead and Anderson (2006) investi- 
gated signatures of stratospheric ozone recovery from MOD 
total ozone for 1979-2005. Their Fig. 3 is similar to Fig. 15 
and showed a similar signature of turnaround in trends over 
the final 10 yr of their ozone record (1996-2005). They con- 
cluded that there was a signature of recovery occurring but it 
may have been influenced by other natural variability involv- 
ing changes in temperature and dynamical forcing. A num- 
ber of other studies have suggested that the early record of 
the turnaround in trends was influenced by natural increase 
caused by change in dynamics, particularly in the Northern 
Hemisphere extra-tropics (e.g., Krizan and Lastovicka, 2005; 
Harris et al., 2008; WMO, 2011). Stolarski and Frith (2006) 
showed a figure (their Figure 7) which is similar to the left 
panel of our Fig. 16 except that they instead plotted monthly 
mean MOD total column ozone and for a shorter time record 
of October 1978-June 2006. Their Fig. 7 also shows an ob- 
vious positive correlation of total ozone with solar UV for 
the Nimbus 7 and Earth Probe periods and with a “level- 
ing off” period of ozone in the final years beginning in the 
late 1990’s despite solar UV decreasing substantially dur- 
ing that time period. An obvious implication of this level- 
ing off period during solar UV decrease is that there may be 
large stratospheric ozone recovery occurring following the 
late 1990’s which overwhelms downward solar forcing and 
dynamical effects. Now with several more years of additional 
ozone measurements than previous studies, stratospheric col- 
umn ozone and total column ozone in Fig. 16 (right panel) 
continue to show increases. Over the next few years solar 
UV forcing will reach maximum and should tend to increase 
global ozone levels further. It will be interesting to see in 
the next few years whether or not global stratospheric ozone 
continues to increase and by how much. 

Among the first papers to evaluate global ozone trends and 
recovery turnaround in models and satellite measurements 
was Jackman et al. (1996) which used a 2-D chemical trans- 
port model. Interestingly, their model predicted a turnaround 
in the mid- 1 990 ’s and recovery of ozone by 2010 similar 
to the amounts present in the mid-1980’s. These features 


are both present in the measurements in the right panel of 
Fig. 16. Similar figures comparing total column ozone from 
models and measurements have been shown recently in the 
2010 WMO report and by Oman et al. (2010, and references 
therein). Most of the recent models predict a turnaround oc- 
curring around year 2000 with largely varying recovery time 
periods. Many of the models predict a longer ozone recovery 
to 1980’s levels by year 2020 or later. 

The Montreal Protocol was an international treaty initiated 
in 1987 to reduce worldwide ozone destroying substances 
such as chlorofluorocarbons and Bromine compounds. As 
a result of this treaty, since 1989 there has been a dra- 
matic global reduction in these substances and a turnaround 
in ozone trends. After accounting for solar UV forcing the 
32-yr record indicates a turnaround in stratospheric ozone 
loss in the mid-1990’s with recent ozone levels comparable 
to amounts present around year 1980. These measurements 
seem to suggest a faster stratospheric ozone recovery occur- 
ring than predicted by many of the recent models. 


6 Summary and conclusions 

The convective-cloud differential (CCD) method is the most 
basic form of the “cloud slicing” technique. This method 
generates gridded measurements of TCO and SCO in the 
tropics and also Pacific mean and zonal mean TCO and SCO 
extending to high latitudes. The CCD method is assessed us- 
ing Aura OMI and MLS ozone measurements for the period 
October 2004 through July 2010. TCO and SCO from the 
CCD method are compared with corresponding OMI/MLS 
residual measurements in which MLS ozone profiles are 
used to calculate SCO. The Aura OMI/MLS data provide 
greatly improved temporal and spatial coverage for evalu- 
ating the CCD measurements as previous validation efforts 
involved sparse data from SAGE and HALOE (for SCO) and 
ozonesondes (for TCO). 

The analyses indicate that the CCD and OMI/MLS grid- 
ded TCO products in the tropics agree to within one DU off- 
set with a mean RMS difference of about 3 DU. RMS dif- 
ferences between these two datasets for inter-annual varia- 
tions in the Pacific are about 0.5-2. 5 DU for either TCO or 
SCO over most latitudes ranging from the tropics to high lat- 
itudes. (RMS differences for TCO are the same as for SCO 
since total column ozone from OMI is used for both resid- 
ual methods.) RMS differences for either zonal mean TCO 
or zonal mean SCO vary from about 0.5-1 DU in the tropics 
to ~ 1.5 DU at mid-to-high latitudes. 

The Aura OMI/MLS ozone data have given us greater con- 
fidence in the CCD measurements of tropospheric and strato- 
spheric ozone derived from previous TOMS measurements. 
This is important for extending the TOMS data with more 
recent OMI data to develop a long record of tropospheric 
and stratospheric ozone for evaluating decadal changes. We 
have developed from TOMS and OMI measurements a long 
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32-yr (1979-2010) dataset of tropospheric and stratospheric 
ozone. The analyses of these time series show that the quasi- 
biennial oscillation (QBO) is the dominant source of inter- 
annual variability of stratospheric ozone. During the Aura 
record the QBO variability in stratospheric ozone was of a 
24-month periodicity and was clearest throughout the South- 
ern Hemisphere extending to high latitudes. Both CCD and 
MLS measurements show QBO-related inter- annual variabil- 
ity of stratospheric ozone in the Southern Hemisphere of 30- 
40 DU. The 32-yr dataset also indicates a QBO signal in tro- 
pospheric ozone with peak-to-peak amplitudes varying from 
about 2 DU up to 7 DU. These results further strengthen evi- 
dence of a QBO signal in tropospheric ozone as suggested in 
previous studies based upon much shorter time record mea- 
surements from SHADOZ ozonesondes and TOMS satellite 
ozone. 

Trend analysis of the 32-yr record dataset suggests a 
turnaround in global stratospheric column ozone loss in the 
mid- 1990’s with current ozone levels comparable to the mid- 
1 980’s. After accounting for solar UV forcing the ozone lev- 
els are now similar to about year 1980. Our analyses show 
that most of the trend turnaround occurred in the Northern 
Hemisphere extending from the subtropics to high latitudes. 
The last 15 yr of record (1996-2010) indicates nearly oppo- 
site positive trends in both hemispheres compared to the first 
15 yr of record (1979-1994). The measurements suggest that 
stratospheric ozone levels comparable to 1980 or earlier may 
be possible within the next several years as the 11 -yr solar 
cycle UV reaches peak production of stratospheric ozone. 

The trends of stratospheric ozone are generally consis- 
tent with the prediction of past and current chemistry-climate 
models, however many of these models predict a turnaround 
occurring around year 2000 and longer ozone recovery to 
mid-1980’s levels by year 2020 or later. Regardless of dis- 
crepancies, our results suggest substantial recovery now oc- 
curring in global stratospheric ozone largely as a result of the 
Montreal Protocol. 
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